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Abstract 

We investigate the avalanche dynamics of the Abehan sandpile 
model on arbitrarily large balls of the expanded cactus graph (the 
Cay ley graph of the free product Z3 * Z2). We follow the approach 
of Dhar and Maiumdar |DM90] to enumerate the number of recur- 
rent configurations. We also propose a method of enumerating all 
the recurrent configurations in which adding a grain to a designated 
origin vertex (far enough away from the boundary vertices) causes 
topplings to occur in a specific cluster (a connected subgraph that is 
the union of cells, or copies of the 3-cycle) within the first wave of 
an avalanche. This filling method lends itself to combinatorial eval- 
uation of the number of recurrent configurations in which a certain 
number of cells topple in the first wave of an avalanche starting at the 
origin, which are amenable to analysis using well-known recurrences 
and corresponding generating functions. Using this, we show that the 
cell-wise first-wave critical exponent of the expanded cactus is | . 
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Science Foundation Grant DMS# 0700811 for research during the summer of 2011. Ad- 
ditionally, the author thanks Tatiana Smirnova-Nagnibeda and Michel Matter for their 
helpful discussions and suggestions, as well as the referees for their review of this paper. 
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1 Introduction 



1.1 Definition of the Abelian sandpile model 

Here we introduce our definition of tlie Abelian sandpile modelj^ 

Let {Gi)'^i be an infinite sequence of finite graphs that is increasing 
(that is, each graph in the sequence is a subgraph of its successor), where 
the sequence exhausts an infinite d-regular transitive graph F (that is, F = 
Ui^i Gi). For example, (Gj)^^ might be a sequence of finite connected graphs 
that each contain a designated origin vertex o in F|^ Hereafter, we will 
consider (Gj)^^ to be an increasing sequence of connected induced finite 
subgraphs of F containing the origin vertex. If we define = Sr(o, = 
{v e V(T)\dr{o,v) < i} (the balls centered at o of radius i according to 
dr, the graph distance function), then the Gj's are an increasing sequence of 
graphs that exhaust F. 

The Abelian sandpile model on each graph Gj in the sequence consists of 
configurations, which are functions x '■ V{Gi) — Z"*" that assign each vertex 
a height that is a positive integer number of grains stacked on that vertex. 
A stable configuration is a configuration where each vertex has height less 
than or equal to d. Transitions from one stable configuration x to another 
are governed by the following rules: 

1. A grain of sand is added to a randomly selected vertex f of G. 

2. If the addition of that grain of sand causes the vertex v to have a height 
exceeding d, then v topples, causing it to lose d grains and causing each 
neighbor of f in G to gain one grain. (There may be fewer than d 
neighbors of v in G, in which case the total number of grains in the 
system will decrease.)]^ 

3. Vertices gaining grains from toppling can also topple themselves if their 

^Good sources for basic information on the Abelian sandpile model found in this section 
are |Dha99j and [MNllj 

^If r is transitive, it does not matter which vertex is designated as the origin. 
Matter and Nagnibeda |MNllj alternatively consider only finite graphs G whose 
vertices can be classified as dissipative; dissipative vertices do not topple. Presumably the 
remaining non-dissipative vertices each have degree d. Given our notion of the Abelian 
sandpile model, we can add the vertices in V{T) — V{G) adjacent to a vertex in G as 
dissipative vertices; the model thus generated is clearly equivalent. 
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height now exceeds d, causing a chain reaction of topphngs 



4. Eventually, the toppling must stop in a configuration where each ver- 
tex's height is less than or equal to d. (This is true since there is at 
least one vertex in G whose degree is less than d; see |Dha99] .) The 
process of toppling vertices iteratively in a configuration to reach a sta- 
ble configuration is known as relaxation. At this point, another grain 
of sand is added to a randomly selected vertex of G. 

If, when we start with a configuration x of G, adding a grain to a vertex 
f of G causes an avalanche, then we can consider relaxing it by toppling a 
sequence of not necessarily distinct vertices Vi, . . . ,Vk- We denote by Xo the 
configuration x with an extra grain at v and Xi the configuration reached from 
Xi-i by toppling t>j. In order for this to be well-defined, we need Xi-ii'Vi) > d 
for all i from 1 to k. We thus define the sequence wi, . . . , f ^ to be a toppling 
sequence if it satisfies the following conditions: 

1. For all integers i from 1 to k, Xi-ii^i) > d, and 

2. The configuration reached after adding a grain to v and toppling vi, . . . ,Vk 
in order is stable. 

Note that, in the sandpile model, unstable configurations never appear 
outside of the process of adding a grain and relaxing. Therefore, from this 
point forward, all configurations are assumed to be stable. 

1.2 The sandpile model as a Markov process 

We can naturally define for each v G V{G) an operator that acts on con- 
figurations by adding a grain to v and relaxing until a valid configuration 
is reached. |Dha99j If we look only at recurrent configurations (described be- 
low), then we can also define inverse operators on configurations x &s 
the unique recurrent configuration x' that is changed to x under a^. The 
operators and their inverses together generate a finite Abelian group. 

If we assign to each v a positive probability then application of a^^s 
successively to recurrent configurations according to those probabilities is a 
Markov process. 

''in the Abelian sandpile model, the order in which vertices topple does not matter. 
This is also true in more generalized sandpile models; see |Dha99) . 
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1.3 Avalanches, recurrent and forbidden configurations, 
and critical exponents 

The series of topplings caused by the addition of a grain to a vertex is called 
an avalanche, and its mass equals the number of vertices that topple at least 
once due to the added grain. 

Within a configuration, a forbidden subconfiguration (FSC) is a set V C 
V with the property that every vertex v & V has a height less than or 
equal to the number of edges incident to v and a vertex of V (see jDM90] ). 
For example, two adjacent vertices with height 1 form an FSC, as do three 
vertices in a chain with heights 1-2-1. [DM90j A configuration is recurrent if 
it almost certainly occurs infinitely often in the Markov process. [ MRZOl] We 
will use the well-known characterization of recurrent configurations. 

Lemma 1 ( |MRZOl] ). A configuration is recurrent if and only if it does not 
have a forbidden subconfiguration. 

In fact, it is possible to equivalently define a configuration x as recurrent 
if it is almost certain that x will appear infinitely many times in the Markov 
process. IMRZOlf l 

In the Abelian sandpile model with sinks (or non-dissipative vertices), 
one can use the burning algorithm to determine whether a configuration x 
is recurrent or contains a forbidden subconfiguration. |Dha99j At the start, 
only the sinks are burned. At each step, the burning algorithm examines 
each vertex v & G and determines whether x(f ) is greater than the number 
of unburnt neighbors — if so, v is burned. This continues until every vertex 
is burned, indicating that x is a recurrent configuration, or until it is impos- 
sible to burn any more vertices (in which case the unburnt vertices form a 
forbidden subconfiguration). |Dha99] The burning algorithm also facilitates 
a bijection between recurrent configurations and spanning forests where each 
component contains exactly one sink. Vertices burned on the first step must 

^Although the Markov process allows sand to be added to any vertex in a given finite 
graph Gi in the sequence of balls of F centered at o, it is valid to only examine what 
happens when a grain is added to o. To examine what happens when a grain is added to 
a different vertex v, a different graph sequence of balls centered at v will need to be used. 
But if r is transitive (as the expanded cactus examined in this paper is) , the results will 
be the same regardless of which vertex is chosen as the origin. Although it is possible to 
consider the avalanche dynamics for every vertex of a fixed finite graph G (which requires 
taking into account boundary effects), this is beyond the scope of this paper. 
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have fewer than d unburnt neighbors and therefore are adjacent to sinks, 
while for a vertex v to burn on the tth step (with t > 1) requires that v not 
be adjacent to a sink and that a neighbor of v burn on the t — 1st step. In a 
recurrent configuration, every vertex burns, and so one can construct a forest 
where each component contains exactly one sink by connecting each vertex 
that burns on the first step to a neighboring sink and each vertex that burns 
on the tth step (t > 1) to a neighbor that burns on the t — 1st step. Although 
there may be more than one possible choice of neighbor, with appropriate 
conventions on which neighbor to select based on the height of the vertex, 
this construction serves as a bijection between recurrent configurations and 
spanning forests where each component contains exactly one sink. |Dha99] 

It is well-known that, in the steady state, each recurrent configuration is 
equally probable. For each k G Z"*", n G Z"*" U {0}, we can calculate both 
the number of recurrent configurations of Gk as well as the number of those 
recurrent configurations in which adding a grain to o causes exactly n distinct 
vertices to topple. Define N^. to be the number of recurrent configurations 
of Gk and Nk^i to be the number of recurrent configurations of Gk in which 
adding a grain to o causes exactly i distinct vertices to topple. Then define 



As k goes to infinity, the values pk{n) should converge for each fixed n. Define 

p{n) = lim Pk{n) 

k—^oo 

Then we can consider the asymptotic behavior of the sequence {p{n)}'^^Q. In 
particular, if there exist constants Ci, C2 > 0, G Z+, and a G M such that 
Gin~°' < p{n) < C2n~°', for all n > A^ where p{n) > 0, then a is the critical 
exponent of the increasing graph sequence {Gi)°Zi exhausting F. 

As we will see, the problem becomes more amenable when we consider 
only those vertices that topple on the first wave of an avalanche. A vertex v 
topples in the first wave of an avalanche caused by adding a grain at if there 
exists a toppling sequence containing v in which the vertices in the sequence 
up to and including the first occurrence of v include exactly once. Define 
N^- to be the number of recurrent configurations of Gk in which adding a 
grain to o causes exactly i distinct vertices to topple in the first wave. Then 
pl{i) and p-^{n) are defined analogously to Pk{i) and p{n). This then allows 
us to state the first-wave critical exponent as a number satisfying that 
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Cin~°' < p^{n) < C2n~°' for all n > where p^{n) > for appropriate 
constants Ci, C2 > 0. 

In this paper, we examine the Abelian sandpile model on arbitrarily large 
but finite subgraphs of the expanded cactus (as named by |FE61] ). which is 
formed by taking the 3-regular Bethe lattice and "decorating" each vertex 
by replacing it with a three- vertex cycle, and the three edges incident to that 
vertex in the Bethe lattice each become incident to a separate vertex in the 
three- vertex cycle. (See Figure flj) 




Figure 1: The infinite expanded cactus with the origin vertex labeled. 

Each copy of the three-vertex cycle in the expanded cactus is designated 
a cell. We will only consider finite expanded cacti that can be generated by 
decorating a finite subgraph of the Bethe lattice (the underlying graph). The 
expanded cactus is the Cayley graph of the group Z3 * Z2 with respect to the 
standard generating set. The maximum height of each vertex is 3. In the 
infinite expanded cactus F, we will take the graph sequence (Gj) defined as 
follows: take the infinite 3-regular Bethe lattice H and fix an origin vertex 
z in H. Gi will be given by the decoration of Bh{z,i), the ball of radius i 
about z in the Bethe latticej^ The cell formed by decorating 2;, present in 

^Matter and Nagnibeda |MNllj examined the critical exponent for cactus graphs in 
the random weak limit. Their hmiting process involves taking larger and larger balls that 
exhaust the graph, but the origin is a randomly selected vertex within each ball. Therefore, 
in the limit, the origin is not a fixed vertex arbitrarily deep within the graph, but rather 
has its origin located according to a probability distribution such that the origin is a finite 
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each Gj, will be designated the origin cell; additionally, one arbitrary vertex 
in that origin cell will be designated the origin vertex. 

Although calculating the distribution of avalanches of specific mass re- 
quires calculating the probability that adding a grain will cause a given num- 
ber of vertices to topple, it is simpler, in the case of the expanded cactus, 
to consider a different definition of avalanche mass: the number of cells in 
which at least one vertex topples (which we will call the cell mass of the 
avalanche). To examine the behavior of the Abelian sandpile model on F, we 
will consider a graph sequence {Gi)°Zi of ^^e expanded cactus F, where Gi is 
the decoration of the ball centered at the origin of the underlying 3-regular 
tree with radius i. Analogously to how Pkin) and p{n) are defined in general, 
we will define ^^(n) and p'^{n) using the number of recurrent configurations 
where exactly n cells topple when a grain is added to the origin vertex in 
place of the number of recurrent configurations where exactly n vertices top- 
ple when a grain is added to the origin vertex. Likewise, we say that a cell 
is toppled on the first wave if any of its vertices is toppled in the first wave; 
then, we define p'^j {n) and p'^^ [n) in a similar fashion using the cells that 
topple in the first wave. 

If there exist constants Ci, C2 > 0, G Z+, and a G M such that Cin^" < 
p^'^in) < 6*2^"", for all n > A^ where p{n) > 0, then a is the cell-wise first- 
wave critical exponent of the increasing graph sequence (Gj)^^ exhausting F. 
If the critical exponent, (respectively, first-wave critical exponent, cell- wise 
first-wave critical exponent) is a for every increasing graph sequence (Gj)^;^ 
exhausting F, then we say that the critical exponent (respectively, first-wave 
critical exponent, cell-wise first-wave critical exponent) of F is a. 

The critical exponent of increasing graph sequences (Gj)^^ exhausting 
the expanded cactus in the Abelian sandpile model is important to mean- 
field theory, which concerns itself with the size of avalanches in critical 
systems. |ZLS95] Dhar and Majumdar |DM90] showed that the critical mean- 
field exponent for infinite binary (3-regular) trees is |. Although it is con- 
jectured that the critical exponents of quasi-isometric graphs should be the 
same, this has not been proven for non-trivial pairs of infinite graphs. Our 
result, below, provides evidence consistent with this conjecture: 

distance away from a vertex with degree less than d with probabihty 1. Thus, the graphs 
studied by Matter and Nagnibeda are one-ended, which is used in their analysis. However, 
the graphs studied in this paper are infinitely-ended, so the approach used in [MNlT] is 
not applicable to our problem. 
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Theorem 1. On the expanded cactus, the cell-wise first-wave critical expo- 
nent is |. 

In percolation theory, there is a similar proposition that the critical ex- 
ponents of percolation on a transitive graph only depend on the large-scale 
structure of the graph |Koz08] and therefore, quasi-isometric graphs will have 
the same critical exponents for percolation. However, there are few calculated 
examples that support this conjecture. 

1.4 Outline of method of Dhar and Majumdar and our 
approach 

Dhar and Majumdar |DM90] determine pi by counting, for each sufficiently 
large subgraph T of the 3-regular tree, the number of recurrent configurations 
(which appears in the denominator) and the number of recurrent configura- 
tions where exactly i vertices topple (which appears in the numerator). Using 
the notions of strongly and weakly allowed configurations on subtrees, Dhar 
and Majumdar first compute the limiting ratio of strongly allowed to weakly 
allowed configurations as the size of a subtree goes to infinity. They then 
derive for every connected induced subgraph {cluster) C with n vertices a 
recurrence relation for the number of recurrent configurations on T in terms 
of the number of strongly allowed configurations of each of the n + 2 subtrees 
Ui, . . . ,Un+2 induced by removing C. In particular, Dhar and Majumdar 
showed that the number of recurrent configurations is 3 ■ 4:"' YYj=i -^siUj) , 
where Ns{Uj) denotes the number of strongly allowed configurations on Uj. 
Then, Dhar and Majumdar show that a configuration x on T is recurrent 
and has the property that adding a grain to the origin causes exactly the n 
vertices in a cluster C to topple if and only if 

1- x{fj) = 3 for every t> G C 

2. Each of the + 2 induced subtrees Ui, . . . , Un+2 has a root vertex whose 
height is 1 or 2 

3. At least one of the n + 2 induced subtrees is strongly allowed 

Dhar and Majumdar show that for each such cluster C of size n, the num- 
ber of recurrent configurations where exactly the vertices in C topple is 
(1 _ 2-('^+2)) n"=f ^s(^i), and this is independent of the shape of C. Dhar 
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and Majumdar then derive the formula for the number of n-vertex clusters 
containing the origin using generating functions, and they showed that the 

3 

number of clusters grows asymptotically to 4"n~2. From this, it follows that 

Pi= (l-2-("+2))^-i 

from which it is clear that the critical exponent of the 3-regular tree is |. 

The remainder of our paper analyzes the cell-wise first-wave critical ex- 
ponent of the expanded cactus using a refinement of the method of Dhar and 
Majumdar. In section 2, we describe how configurations can be enumerated 
by dividing the finite graph into decorated rooted subtrees and defining and 
considering weakly and strongly allowed radicals (configurations restricted 
to decorated rooted subtrees) on those decorated subtrees. Based on this 
method, we can express the number of recurrent configurations on a finite 
graph with regards to a cluster (defined below) of cells about the origin in 
terms of the number of strongly allowed radicals on the decorated subtrees 
induced by removing the cluster. Section 4 outlines how the number of recur- 
rent configurations on a rooted expanded cactus where cells in a given cluster 
topple and no others can be enumerated in terms strongly and weakly allowed 
radicals induced by removing the cluster, using the concept of filling rules 
to determine what configurations are allowed within the cluster. Section 5 
examines how we can adapt the filling rules to enumerate configurations on 
the expanded cactus where given cells topple on the first wave, and it ex- 
amines issues with multiple-wave avalanches on the expanded cactus, issues 
arising from our graph's being a cactus and not a tree as in |DM90j . In 
section 6, we use asymptotic enumeration methods to express the number 
of configurations where n cells topple in terms of the number of strongly 
allowed radicals induced by any cluster of size n. Our method of counting 
overcounts in that it includes configurations that are not recurrent — section 
7 shows that only a small number of positions are erroneously counted and 
thus that the asymptotics are not affected. Finally, section 8 provides open 
questions and a framework for generalizing the method used in this paper to 
arbitrary decorations of infinite trees by finite transitive connected graphs. 
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2 Counting method — weakly and strongly al- 
lowed radicals 



Throughout this paper, our method of enumerating recurrent configurations 
on the finite expanded cactus G with certain properties is to take a cluster (or 
connected subgraph formed by the union of cells) C with n cells, containing a 
predetermined cell called the origin, one of whose vertices is also designated 
the origin vertex. The subgraph induced hy V \ C has n + 2 connected 
components, which are decorated rooted subtrees. Each decorated rooted 
subtree has a cell immediately adjacent to the cluster, which we will call the 
root cell, and a root vertex within the root cell that is adjacent to a vertex 
within the cluster. 

Given a configuration x and a decorated rooted subtree f/, the restriction 
xljy is called a radical. The radical's root cell and vertex are the same as the 
decorated rooted subtree on which it is based. Considering the radical as not 
attached to anything, the configuration may be a recurrent configuration or 
have an FSC If the radical does not have a forbidden subconfiguration, then 
we will create a dichotomy, analogous to the one introduced in |DM90j . of 
strongly allowed and weakly allowed radicals. A radical that does not contain 
an FSC is strongly allowed if, when attached to a vertex of height 1, does not 
result in an FSC, and is weakly allowed otherwise. Also, by definition, if U is 
empty, then the null configuration on f/ is a strongly-allowed radical. Given a 
decorated rooted subtree U , we denote the number of weakly allowed radicals 
on U as Nyu{U), the number of strongly allowed radicals on U as Ns{U), and 
the ratio ^^f^ as Xu. 

In our enumerations, our goal is to express the number of recurrent config- 
urations, as well as the number of recurrent configurations wherein topplings 
occur in exactly the cells of a given cluster C of n cells, in terms of the 
product 

n+2 
i=l 

where each Ui is a decorated rooted subtree induced by C. We first turn to 
counting the number of recurrent configurations. 
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3 Enumerating the recurrent configurations 



Consider the origin cell in the expanded cactus and a configuration of both 
the cell and the three induced radicals. Whether the configuration is recur- 
rent or not can be determined entirely by the heights of the vertices in the 
origin cell, as well as whether each radical is strongly or weakly allowed. 

In Table [T} we examine the 16 possible cell configurations and determine 
which combinations of strongly and weakly allowed radicals result in an re- 
current configuration. Here, S denotes a strongly allowed radical, W a weakly 
allowed radical. The order of the letters in each combination indicates what 
type of radical is attached where, corresponding to the order of the vertices 
in the given cell. For example, the entry WSS in the 2-2-3 row indicates that 
if a weakly allowed radical is attached to the first 2 and strongly allowed 
radicals to the second 2 and to the 3, the resulting configuration is recurrent. 



Cell recurrent radical combinations 



1-2-3 


sss 








1-3-2 


sss 








2-1-3 


sss 








2-3-1 


sss 








3-1-2 


sss 








3-2-1 


sss 








2-2-3 


sss. 


WSS, 


sws 




2-3-2 


sss. 


WSS, 


ssw 




3-2-2 


sss. 


sws. 


ssw 




1-3-3 


sss. 


sws. 


ssw 




3-1-3 


sss. 


WSS, 


ssw 




3-3-1 


sss. 


WSS, 


sws 




2-3-3 


sss. 


WSS, 


sws. 


ssw, wws, wsw 


3-2-3 


sss. 


WSS, 


sws. 


ssw, wws, sww 


3-3-2 


sss. 


WSS, 


sws. 


ssw, wsw, sww 


3-3-3 


sss. 


WSS, 


sws. 


ssw, wws, wsw, sww 



Table 1: The combinations of allowed radicals for each allowed origin cell 
configuration that result in a recurrent configuration. 

We also must consider how to compute the number of strongly and weakly 
allowed radicals on a given decorated rooted subtree. This can be calculated 
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inductively by considering the 16 possible configurations of the root cell of a 
radical, as well as whether the two induced subradicals are weakly or strongly 
allowed. Each combination of a root cell configuration and weakly or strongly 
allowed subradicals produces either a radical with a forbidden subconfigura- 
tion, a weakly allowed whole radical, or a strongly allowed whole radical. 

In Table [2| we illustrate the combinations of root cells and weakly and 
strongly allowed radicals that produce weakly and strongly allowed whole 
radicals. The root vertex is the last number listed in the cell; the two letter 
codes indicate whether a weakly allowed (W) or strongly allowed (S) radi- 
cal is being added to the vertex indicated by the first and second number, 
respectively. 



Cell 


Weakly allowed whole 


Strongly allowed whole 


1-2-3 


SS 




1-3-2 


SS 




2-1-3 


SS 




2-3-1 


SS 




3-1-2 


SS 




3-2-1 


SS 




2-2-3 


SS, WS, SW 




2-3-2 


ws 


SS 


3-2-2 


SW 


SS 


1-3-3 


SW 


SS 


3-1-3 


ws 


SS 


3-3-1 


SS, ws, SW 




2-3-3 


SW, WW 


SS, ws 


3-2-3 


ws, WW 


SS, SW 


3-3-2 




SS, ws, SW 


3-3-3 


WW 


SS, ws, SW 



Table 2: The combinations of allowed radicals for each allowed root cell 
configuration that result in strongly allowed or weakly allowed whole radicals. 

From Table [2| we can calculate the limiting ratio of weakly allowed to 
strongly allowed radicals on arbitrarily large decorated rooted subtrees. To 
do this, we first tabulate the number of instances of WW, WS, SW, and SS 
in the weakly allowed and the strongly allowed columns: 



12 





WW 


WS 


sw ss 


Weakly allowed whole 


3 


5 


5 8 


Strongly allowed whole 





3 


3 8 



Table 3: The number of root cell configurations for which the given combi- 
nations of allowed subradicals produce a weakly allowed or strongly allowed 
whole radical. 

Suppose [/ is a decorated rooted subtree such that when the root cell 
is removed, two smaller decorated rooted subtrees, Ui and U2 are induced. 
From Table [3} we know 

N^{U) = 3N^{Ui)N^{U2) + 5N^{Ui)Ns{U2) 
+ 5NsiU,)NUU2) + 8N,iU^)N,iU2) 
= {3xu,xu, + 5xu, + 5xu, + 8)Ns{Ui)Ns{U2) 

and 

N,{U) = 3N^{Ui)Ns{U2) + 3Ns{Ui)N^{U2) + 8Ns{Ui)Ns{U2) 
= {3xu, + 3xu, + 8)Ns{Ui)Ns{U2) 

Therefore, 

_ 3xuiXu2 + 5xui + 5xu2 + 8 
3xui + 3xu2 + 8 

Consider now the decorated rooted subtrees generated as follows: Bq as a 
single cell, i?„ as a decorated rooted subtree whose child subtrees are -B„_i. 
Let Xn = j^*"^^^"! for each nonnegative integer n. Then 

Xo = l 

and 

^ (8 + 10x„„i + 3x^_i)Ar,(T„„i)2 ^ 8 + 10x„.-i + 3xl_^ 
(8 + 6x„_i)iV,(T„_i)2 8 + 6x„_i 

This simplifies to 

1 
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so X, 



n 



— 2 — 2 and it is clear that 



lim x< 



n 



2 



Further, we can compute 



dxu 
dxu^ 



as follows: 



dxui 



dxu 



(3x2 + 5) (3x1 + 3x2 + 8) — 3(3,xi,X2 + 5xi + 5x2 + 8) 
(3x1 + 3x2 + 8)2 



9x1 + 24x2 + 16 
(3xi + 3x2 + 8)2 



This means that xu increases as a function of xu^ and, by symmetry, Xu2- 
Thus, if ?7 is a decorated rooted subtree that contains a subtree B^, then xu > 
Xn — 2 — 2~". Further, since xc/(2,2) = 2, it is trivial to show inductively 
that Xu < 2 for every finite decorated rooted subtree U . Together, these 
results show that we can guarantee that a decorated rooted subtree U will 
have xy arbitrarily close to 2 simply by ensuring that IJ is deep enough. 

Further, we have the following result that relieves us from having to 
consider which specific graph sequence we are using: 

Theorem 2. Let [Gn)'^^^ '^^ increasing graph sequence that exhausts T, 
the infinite expanded cactus, where each Gn is a connected graph formed by 
the union of cells, including the origin. Then, for every e > 0, and for every 
positive integer k, there exists a positive integer N such that for every cluster 
C of at most k cells in containing the origin cell and every decorated 
rooted subtree U induced by G, it holds that 2 > xu > 2 — e. 

Now let us suppose that C is a cluster of n cells about the origin, in- 
ducing decorated rooted subtrees Ui, . . . , f/„+2- We would like to express the 
number of recurrent configurations on the entire graph in terms oi Ns{Ui) ■ 
Ns{U2) ■ ■ ■ Ns{Un+2)- The following analysis works with clusters of any shape, 
but for simplicity, we will assume that the cluster is a chain of n cells, num- 
bered consecutively from 1 to n. Label the two induced decorated rooted 
subtrees adjacent to cell 1 Ui and U2, the one decorated rooted subtree ad- 
jacent to cell k (for 1 < k < n) Uk+i, and the two decorated rooted subtrees 
adjacent to cell n Un+i and Un+2- Also, for 1 < A; < n, define to be the 
decorated subtree whose root is cell k and which does not contain cell A; -|- 1. 

With these labellings, wc can compute the total number of recurrent 
configurations on the graph by dividing the graph at cell n and using Table 
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[H The result is 

A^recurrent = [16 + 8Xt„_i + 8Xu„+i + 8Xu„+2 
Ns{Tn-l)Ns{Un-,l)Ns{Un+2) 

Further, we have the recurrence, for 2 < k < n, 

NsiW = [8 + 3xu,^, + 3xT,_JiV,(f/,_,0iV.(Tfc-i) (1) 

and 

iV,(Ti) = [8 + 3xu, + 3xu2]NsiU,)NsiU2) (2) 

If e > is given and G is a finite graph satisfying the conclusion of 
Theorem |2j then we have 

iVrecurrent > [100 - 60e + 96^] iV, (r„_i)iV, ([/„+i)iV, (f/^+s) 

and iterative expansion gives 

n+2 

A^recurrent > (20 - 6e)"-^(100 - 60e + 96^) JJ Ns{U,) 

i=l 

Also, since xu < 2 for every finite decorated rooted subtree U, we have 

n+2 

A^recurre„t< 100-20"-lJ]iV,(f/,) 

i=l 

Therefore, as e approaches 0, we have 

n+2 

A^recurrent = 5-20" J]iV,([/,) 

i=l 

4 Enumerating the recurrent configurations 
on a decorated rooted subtree in which cells 
of a given cluster will topple using filling 
rules 

In this section, we propose filling rules that, given a decorated rooted subtree 
and a cluster C about the origin, describe what combinations of configura- 
tions of the cells of C and radicals induced by C will result in an recurrent 
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configuration wherein adding a grain to the origin causes topphngs to occur 
within at least one vertex in each cell of C and no vertices outside the cells 
of C. The rules allow for independent assignment of configurations to each 
cell of C, and requirements for the radicals only depend on the cell to which 
the underlying decorated rooted subtree is attached. Therefore, the number 
of recurrent configurations on a decorated rooted subtree where a given clus- 
ter of cells topples can be enumerated by simple use of the multiplication 
principle of counting. 

4.1 Terminology 

To explain the filling rules clearly, we need to define some additional terms. 
Within each cell, the origin-facing vertex is the vertex closest to the origin 
vertex. Each cell c within C is classified as internal, medial, or terminal, 
as follows. If each non-origin-facing vertex of c is attached to a cell in the 
cluster, then c is an internal cell. If exactly one non-origin-facing vertex of c 
is attached to a cell in the cluster, then c is a medial cell. If no non-origin- 
facing vertices in c are attached to cells in the cluster, then c is a terminal 
cell. Figure [2] illustrates an example cluster and the internal, medial, and 
terminal cells within it. 




Figure 2: An example cluster, with cells labeled as internal (I), medial (M), 
or terminal (T). 

Define a radical to be a stopper if it is allowed and its root vertex has 
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height 1 or 2. Also, we define the empty radical to be a stopper. For radicals 
that are sufficiently deep, the following hold: 

1. The number of stoppers equals the number of strongly allowed radicals 

2. The number of strongly allowed stoppers is ^ of the number of strongly 
allowed radicals. 

To see this, refer to Table [2j Since weakly allowed subradicals are twice as 
prevalent as strongly allowed subradicals, weight each SS as 1, each SW and 
WS as 2, and each WW as 4. Then calculating the total weights for strongly 
allowed whole radicals, stoppers, and strongly allowed stoppers, gives the 
preceding results. 

4.2 The rules 

The recurrent configurations in which topplings occur in each cell in C and 
no others are precisely characterized by the following: 

1. The cells in C have configurations specified by the filling rules, accord- 
ing to their being internal, medial, or terminal. 

2. The radicals are all allowed, and they satisfy the specific requirements 
for being strongly allowed or a stopper that are determined by the cell 
to which the underlying decorated rooted subtree is attached. 

4.2.1 Internal cells 

The possible configurations for internal cells are 

1. 3-3-3 

2. 3-3-2, with the vertex closest to the origin having height 3 (this item 
encompasses two cell configurations by symmetry) 

^This holds regardless of the shape of the radical, so long as it is a decorated rooted 
subtree, since there is a bijection between stoppers and strongly allowed radicals. 
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4.2.2 Medial cells 



The possible configurations for medial cells, as well as the corresponding 
restrictions on the attached radical, are 

1. 3-3-3 — radical must be a stopper 

2. 3-3-1, with the vertex adjoining the induced decorated rooted subtree 
having height 1 — radical must be strongly allowed 

3. 3-3-2, with the vertex adjoining the induced decorated rooted subtree 
having height 2 — radical must be a stopper 

4. 3-3-2, with the vertex adjacent to another cell in C but not the closest 
vertex to the origin having height 2 — radical must be a stopper 



4.2.3 Terminal cells 

Finally, we turn to terminal cells, which have these possible configurations 
that can be filled into the cell (and the corresponding radical restrictions): 

1. 3-3-3 — both radicals must be stoppers 

2. 3-3-2, with the vertex closest to the origin having height 3 — both rad- 
icals must be stoppers (this item encompasses two cell configurations) 

3. 3-3-1, with the vertex closest to the origin having height 3 — the radical 
on the decorated rooted subtree adjacent to the 3 must be a stopper; 
the other radical must be strongly allowed (this item encompasses two 
cell configurations) 

4. 3-2-2, with the vertex closest to the origin having height 3 — at least 
one of the two radicals must be strongly allowed; the other must be 
allowed 

5. 3-2-1, with the vertex closest to the origin having height 3 — both rad- 
icals must be strongly allowed (this item encompasses two cell config- 
urations) 
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4.3 Correctness of filling rules for decorated rooted 
subtrees 

In order to use the filling rules to count the number of recurrent configurations 
wherein the cells in cluster C topple and no others, we prove the following 
result. 

Theorem 3. Let x be a configuration on G, a finite decorated rooted subtree, 
and let C be a cluster of G. Then x is a recurrent configuration and adding 
a grain to the origin vertex causes exactly the cells in C to topple if and only 
if X can be generated by applying the filling rules to the cluster C in G. 

To prove this, we will use the following propositions. (Unless otherwise 
stated, G is a connected finite subgraph with origin of an infinite connected 
(i-regular graph F.) 

Lemma 2. Let v E G be given such that x{v) = d. Let vi,V2, ■ ■ ■ ,Vk be a 
toppling sequence for the avalanche caused by adding a grain to v. Then, if 
Vi,V2, . . . ,Vk are not all distinct, then the first repeated vertex in the sequence 
is V. 

Proof. Take t to be the smallest positive integer such that Vf = Vf for some 
t' < t. For each s from 1 to k, let Xs denote the not necessarily stable 
configuration obtained by starting with adding a grain to u, and relaxing 
only the vertices vi, . . . ,Vs in that order. Now, we have that Xt-i{vt) > d 
since Vf is the tth toppling vertex in the sequence. But we also know that 

t-i 

Xt-i{vt) = xi'^t) + ^v) + ^ l{vs r^vt)-d 

However, the sum l{vs ~ ft) is at most d, since each of the d neighbors 
of Vt topples at most once prior to the second toppling of Vt- Thus, in order 
for xt-i{yt) > d to hold, we must have Vt — v. □ 

Lemma 3. Let v & G be given such that x(f ) = d. Then v topples more 
than once in the avalanche caused by adding a grain to v if and only if v 
has d neighbors in G and each of those neighbors topples. In particular, if 
daiy) < d, then no vertex topples more than once in the avalanche caused by 
adding a grain to v. 
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Proof. Let Vi,...,Vk be a toppling sequence for the avalanche caused by 
adding a grain to x 

{=^) Suppose V topples more than once. Let t be the positive integer 
such that Vt = V and such that there exists exactly one t' < t such that 
Vt' — V. Then, for each s from 1 to k, let Xs denote the not necessarily stable 
configuration obtained by starting with adding a grain to u, and relaxing 
only the vertices vi, . . . ,Vs in that order. We have 

f-i 

Xt-i{v) = 0? + 1 + ^ l{vs r^vt)-d 

s=l 

But since v topples at step t of relaxation, we have 

Xt-i{v) >d+l 

This can only happen when X]!=i l('^s ~ '^t) maximum value, d. When 
Yfs^i ~ ''^t) = d, we have that v has d neighbors in G and each of those 
neighbors topples, as desired. 

(«^=) Suppose that each neighbors of v topples more than once. Choose 
t to be the smallest positive integer such that for each neighbor u of v, there 
is a positive integer s <t such that Vs — u. Then, by the converse, v has not 
toppled twice, so it has toppled exactly once by step t. We have 

t-i 

Xtiv)^d + l + J2H'",r.Vt)-d 

s=l 

and by hypothesis, we know that Yll^i ~ '"t) — d, so 

Xt{v) ^d + 1 

and the relaxation process requires that v topple again. 
The second statement of the lemma is immediate. 

□ 

In particular, given a radical of the expanded cactus, adding a grain to 
the root vertex causes each vertex to topple at most once. 

Lemma 4. Let v & G be given such that = d. Let vi, . . . ,Vk be a 
toppling sequence for the avalanche caused by adding a grain to v. For each 
t — 1,2, ... ,k, define St — {vs\l < s <t}. The subgraph of G induced by St 
is connected and contains v. 
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Proof. Let t, 1 < t < A; be given. Clearly, v G S. Let H be the subgraph of 
G induced by St, and suppose for a contradiction that H has more than one 
component. Let C be a component of H that does not contain v. Choose 
the smallest positive integer s such that Vg e St- By this choice of s, we have 
the no neighbor of Vg has toppled by the sth step of relaxation. Therefore, 
Xs-i(fs) = x{vs) < d, so we are not allowed to relax Vs at the sth step, which 
gives us a contradiction. Hence, H is connected. □ 

The following is an immediate consequence. 

Lemma 5. Let G be a finite decorated rooted subtree with root vertex a having 
height 3. Then the cells that topple in the avalanche caused by adding a grain 
to o form a cluster. 

Lemma 6. Let G be a finite decorated rooted subtree with root vertex o with 
configuration x such that x(o) = 3, and let c be cell with origin-facing vertex 
having height less than 3 (in particular, c is not the origin cell). Then c does 
not topple in the avalanche caused by adding a vertex to a. 

Proof. Suppose that the avalanche caused by adding a vertex to a has a 
toppling sequence Vi,...,Vk. Let v be the origin-facing vertex of c. We 
claim that v does not topple. Suppose for a contradiction that Vt = v for a 
(necessarily unique) positive integer t. We note that w is a cut-vertex of G, 
and only one neighbor m of w in G is contained in the connected component 
oiG\{v} containing a. Thus, by Lemma|4| we have that St = {fs|l < s < t] 
induces a connected subgraph of G. Therefore, the only neighbor of v that 
could topple before v does is u. But we then have Xt~i{v) = x{v) + 1 < d + 1, 
so V cannot be the tth vertex to topple. This is a contradiction. So v does 
not topple. The other two vertices in c are not in the same component of 
G \ {v} as o, and so they cannot topple either. Hence, c does not topple. □ 

The following transfer rule allows us to manipulate toppling sequences in 
a way that makes the proof of the main theorem easier. 

Lemma 7. Let G with origin o and configuration x be given such that x{o) = 
d. Suppose that vi, . . . ,Vk is a toppling sequence where no vertex appears more 
than once. Ifl<i<j<k and for all m such that i < m < j , we have 
Vj 7^ f m; then the sequence f i, . . . , Vj, Vi, . . . , f j-i, ^^j+i, ■ ■ ■ ,Vk is a toppling 
sequence. 



21 



Proof. We know that each vertex topples at most once, so each vertex appears 
at most once in the topphng sequence. In the original sequence, the number 
of grains on vj at the time it topples is equal to + St=i ~ ^i)- 

Because none of the vertices are adjacent to Vj, we have that 

the number of grains on vj before Vi topples is also x('^i) + J2l=i ^i'^t ~ Vj), 
which is the same as the number of grains on vj at the time Vj topples. 
Consequently, moving vj before Vi will still permit us to legally topple Vj at 
its new position in the sequence. The other vertices Vi, . . . ,Vj-i will have 
the same number of grains when they topple in the new sequence as they 
did in the old sequence. Finally, the same vertices topple in the same order, 
so after the topplings in each sequence, the configurations are the same. 
Therefore, the new sequence gives a stable configuration after the topplings 
are finished. Thus, the new sequence, vi, . . . ,Vj,Vi, . . . , f j-i, Vj+i, . . . ,Vk, is 
a toppling sequence. □ 

Lemma 8. Let G be a finite decorated rooted subtree with origin vertex a 
with configuration x such that x(o) = 3. Let c be a cell in G, and let v be 
the origin-facing vertex of c. Let e be the edge connected v to a vertex not 
in c, and define U to be the decorated rooted subtree that is the connected 
component of G \ {e} containing v. If v topples in the avalanche caused 
by adding a grain to a, then, for every vertex u G V{U), u topples in the 
avalanche in G caused by adding a grain to a if and only if u topples in the 
avalanche in U caused by adding a grain to v. 

Proof. Let x be a configuration on G, and let vi, . . . ,Vk be a toppling se- 
quence. By hypothesis, v is in this sequence as vertex ve. First, we claim 
that there is a permutation of this sequence such that all the vertices after v 
in the sequence are in V{U). We use Lemma [Tjiteratively as follows: take the 
smallest index i such that i > i and Vi ^ V{U). By Lemma |4| we have that 
the neighbor of v that is not in U appears before v in the toppling sequence. 
Thus, since each vertex topples at most once, we have Vi 7^ v. Moreover, 
Vi ^ w for any w G ViU^-, including . . . , t>j_i. Thus, we can apply 
Lemma [7] to bring Vi directly before vi. After completing this process, we 
have a toppling sequence of the form w\, . . . , Wm-, v,u\, . . . , Um', where each 
Wi ^ V{U) and each Ui G V{U). 

Now, we want to show that v,ui, . . . , Um' is a toppling sequence for the 
avalanche caused by adding a grain to the configuration at f in ?7. For 
each i, 1 < i < m', we have that Wj's neighbors are in V{U). Therefore, Uj's 
height just before it is toppled is the same in both wi, . . . ,v,ui, . . . , Um' and 
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f , Ml, . . . , Um'- Moreover, after toppling f , mi, . . . , m^', the resulting configu- 
ration {x\u)m'+i niust be stable, for if this were not the case, we would be able 
to topple some vertex z E U, which we could have also toppled at the end of 
the sequence Wi, . . . , Wm, v,Ui, . . . , Um'- Thus, the sequence v,Ui, . . . , Um' is a 
toppling sequence for the avalanche caused by adding a grain to configuration 
xlc/ at f in f/, and the lemma follows. □ 

Lemma 9. Let x be a configuration on the decorated rooted subtree G with 
origin a. Suppose x(o) = 3, and let vi, . . . ,Vk be a toppling sequence for the 
avalanche caused by adding a grain to a. If x has a forbidden subconfiguration 
S, then no vertex in the sequence Vi, . . . ,Vk is in S . 

Proof. To show this, it suffices to show that the vertices Vi, . . . ,Vk can be 
sequentially burned according to the burning algorithm. Indeed, vi, which 
must be o, can be burned because = 3 and d{o) = 2. Now, for any 
positive integer i from 2 to n, the number of unburnt neighbors of Vi remaining 
after we have burned fi, . . . ,fj-i is at most 3 — ^{vj ~ fj). But we 

know that xi^^i) + X]j=i ~ "^'i) > 3. This implies that xi^i) > 3 — 
~ ^i)' which is more than the number of unburnt neighbors of Vi. 
Consequently, we can burn Vi at that stage. By the burning algorithm, no 
forbidden subconfiguration 5* can contain any of the vertices Vi, . . . ,Vk □ 

Proof of Theorem^ By strong induction on the number of cells in 

C . In the base case, suppose that C is a one-cell cluster. Then the origin 
cell is a terminal cell. Note that every allowed cell configuration with origin- 
facing vertex having height 3 is allowed to be filled into a terminal cell. We 
therefore check all five cases for the cell configuration (up to permuting the 
two non-origin heights): 

1. 3-3-3 and 3-3-2. The two non-origin vertices topple in the avalanche, 
and so the attached radicals must either be empty or have a root vertex 
with height less than 3. In either case, the radicals are stoppers, as 
prescribed by the filling rules. 

2. 3-3-1. The non-origin vertex with height 3 topples in the avalanche, 
and so the attached radical must either be empty or have a root vertex 
with height less than 3, making such radical a stopper. Moreover, if 
the radical attached to the 1 is not strongly allowed, then it forms a 
forbidden subconfiguration when attached to the 1. Thus, the radical 
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attached to the 1 must be strongly allowed. These are precisely the 
conditions imposed by the filling rules. 

3. 3-2-2. We need to show that at least one of the attached radicals 
must be strongly allowed. For if both are weakly allowed, then the two 
radicals, when joined to the two vertices with height 2, have a forbidden 
subconfiguration. 

4. 3-2-1. We must show that both radicals must be strongly allowed. If 
the radical attached to the 2 is weakly allowed, then that radical, union 
the two non-origin vertices of the origin cell, contains a forbidden sub- 
configuration. If the radical attached to the 1 is weakly allowed, then 
that radical, when attached to the 1, has a forbidden subconfiguration. 

Thus, the base case is complete. 

Now for the induction step, we let C be a cluster with n > 2 cells, and 
we assume that {=^) holds for any cluster with fewer than n cells. Note 
that the origin cell of C must be internal or medial. We handle both cases 
in turn. 

If the origin cell is medial, then exactly one non-origin vertex in the origin 
cell is attached to another cell in C. Call that vertex a, and call the other 
non-origin vertex in the origin cell b. Denote by c the vertex not in the 
origin cell that is adjacent to a. Let U denote the decorated rooted subtree 
consisting of C minus the origin cell, rooted at vertex c. Since c topples, a 
must topple. The only four cell configurations for the origin cell in which 
a topples are 3-3-3, 3-3-2, 3-3-1, and 3-2-3 (where each triple of numbers 
denotes the heights of o, a, and b, respectively). We have four cases to check 
to ensure that the radical attached to b satisfies the filling rules: 

1. 3-3-3, 3-3-2, and 3-2-3. Since b topples, the origin of the radical 
attached to b must have height less then 3 (or else it would topple, 
too). Thus, that radical must be a stopper. 

2. 3-3-1. If the radical attached to b is weakly allowed, then we have a 
forbidden subconfiguration when that radical is attached to vertex b. 
Thus, the radical attached to b must be strongly allowed. 

This means that the origin cell and the attached radical must satisfy the fill- 
ing rules. Moreover, by Lemma [8} when a grain is added to the configuration 
x\u at c, the cells that topple are precisely the cells in the cluster [/ n C of 



24 



n — 1 cells. Thus, by the induction hypothesis, those cells and their attached 
radicals must satisfy the filling rules. Therefore, the entire graph G satisfies 
the filling rules with respect to cluster C. 

Now suppose that the origin cell is internal. Then we must have all three 
vertices in the origin cell topple. This means that the origin cell must be 
either 3-3-3 or 3-3-2 (with the 2 in a non-origin vertex). Removing the origin 
cell induces two subradicals Ui and f/2, and again by Lemma [8} the cells that 
topple in f/i and U2 are precisely C fl t/i and C fl t/2, which are both clusters 
of less than n cells. Those cells and their attached radicals must satisfy the 
filling rules, and so the whole cluster satisfies the filling rules. This completes 
the induction and the first half of the proof. 

(<^=) For the second half of the proof, suppose that G has a configuration 
that satisfies the filling rules with respect to a cluster C in G. By inspection, 
it is clear that all of the cells in C topple. It remains to show that G does 
not contain a forbidden sub configuration. Suppose for a contradiction that 
S C V{G) is a forbidden sub configuration. Then, we can take a connected 
component of S", and that will also be a forbidden subconfiguration. Note 
that S cannot be a subset of a radical induced by C, as each radical must be 
allowed. Also, by Lemma [9} it follows that the vertices in 5* are confined to 
the nontoppling vertices in a single cell c in the cluster and one or both of 
its attached radicals. We examine each case of c: 

1. c is a 3-3-3, 3-3-2, or 3-2-3 medial cell. This is impossible, as all 
of the vertices of c topple. Thus, S* is a subset of the attached radical, 
which is impossible. 

2. c is a 3-3-1 medial cell. This implies that the attached radical is 
weakly allowed, which is barred by the filling rules. 

3. c is a 3-3-3, 3-3-2, or 3-2-3 terminal cell. This is impossible, as all 
of the vertices of c topple. Thus, S* is a subset of the attached radical, 
which is impossible. 

4. c is a 3-3-1 or 3-1-3 terminal cell. Then S must consist of the 
nontoppling vertex with height 1 plus some vertices in the attached 
radical. Thus, the attached radical is weakly allowed, when it is, in 
fact, required to be strongly allowed by the filling rules. 

5. c is a 3-2-2 terminal cell. Then S must contain both nontoppling 
vertices with height 2. Let Ui and U2 be the attached radicals. Then 
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S n V{Ui) and S fl V{U2) have the property that every vertex except 
the root has at least as many neighbors in S as its height. Therefore, 
Ui and U2 would each contain a forbidden subconfiguration if its root 
were attached to a vertex with height 1. Thus, Ui and U2 are both 
weakly allowed, which is not permitted by the filling rules. 

6. c is a 3-2-1 or 3-1-2 terminal cell. We have two possibilities. If the 
nontoppling vertex with height 2 is not in S, then S contains the vertex 
with height 1 and part of its attached radical, making that attached 
radical weakly allowed. But the filling rules do not permit this. If 
the nontoppling vertex with height 2 is in S, then S must contain the 
vertex with height 1 also. Then, the attached radical U to the vertex 
with height 2 has a subset V{U) n S satisfying that every vertex in 
this subset, except for the root, has at least as many neighbors in the 
subset as its height. Therefore, if a vertex with height 1 is attached 
to U, this subset, together with the newly- attached 1, would form a 
forbidden subconfiguration. Therefore, U is weakly allowed. Again, 
though, the filling rules do not permit this. 

Therefore, we cannot find any such forbidden subconfiguration S. Thus, as 
desired, x is recurrent and adding a grain to the origin causes the cells in C 
to topple and no others. This completes the proof. □ 

4.4 Overall impact of filling rules 

The requirements for each cell configuration are independent, and the re- 
quirements for each radical depend only on the cell to which the underlying 
decorated rooted subtree is attached. With the exception of a 3-2-2 termi- 
nal filling, each assignment of a configuration to each cell in C accounts for 
j--j-n+2 j^f^i^jj.-^ configurations, since each radical Ui must be strongly allowed or 
a stopper, either of which is a restriction that allows for Ns{Ui) configurations 
on the radical. 

However, a 3-2-2 terminal filling has a different effect on the number 
of allowed radicals on the two attached decorated rooted subtrees. Since 
neither radical has to be a stopper, but only one has to be strongly allowed, 
this allows for three possibilities: both radicals are strongly allowed, the 
first is weakly allowed and the second strongly allowed, or vice versa. Since 
weakly allowed radicals are twice as prevalent as strongly allowed radicals in 
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the limit, this means that there are 5 times as many configurations for the 
two attached radicals as there would be for any other terminal filling. 

Therefore, we can state that if we have a configuration of the cells of C 
where there are k terminal cells in C with the 3-2-2 configuration, then the 
number of overall configurations accounted for by the given configuration of 
Cis 

n+2 

1=1 

An easy and appropriate way to account for this fact is to treat the 
possibility of a 3-2-2 terminal filling as if it were five separate fillings. If 
we count a 3-2-2 terminal fillings as 5 effective fillings, and all other fillings 
(internal, medial, and terminal) as 1 effective fillings, then we have 

1. 3 effective fillings for each internal cell 

2. 4 effective fillings for each medial cell 

3. 12 effective fillings for each terminal cell 

and the total number of recurrent configurations in which topplings occur in 
each cell in C and no other cells is 

n+2 

fe=i 

where i,m,t are the number of internal, medial, and terminal cells of C. 

5 Extending the filling rules to finite copies 
of the expanded cactus, generally 

In the previous section, we showed that the filling rules enumerate all recur- 
rent configurations on a decorated rooted subtree in which adding a grain to 
the origin causes exactly n cells to topple. However, our main interest is to 
consider all recurrent configurations on a sufficiently large finite copy of the 
expanded cactus where adding a grain to the origin causes exactly n cells to 
topple. Although this extension may seem like it can be easily solved, the 
extension to the expanded cactus presents complications that make enumer- 
ating the configurations much more difficult. 



27 



5.1 The liberty rule 



Let G be a finite expanded cactus with origin o. If C is a cluster about tlie 
origin tliat does not contain tlie cell directly opposite o, then the height of 
the neighbor of o not in the origin cell is less than 3. Thus, no cell topples 
more than once, and the filling rules can be usec^to determine which stable 
configurations have the property that adding a grain to o causes precisely 
the cells in C to topple. 

The issue is that positions resulting from applying the filling rules no 
longer have to be recurrent. To fix this, we impose an additional requirement 
on our fillings. In a configuration x define a liberty to be a radical induced 
by C that is attached to a vertex in C that is connected to o by a path 
of vertices of height 3. Note that the number of liberties depends on the 
configuration. We then require that at least one liberty be strongly allowed]^ 
This is justified by the following lemma. 

Lemma 10. Let G be a finite expanded cactus with origin o, and let C be 
a cluster of G about the origin that does not contain the cell directly oppo- 
site 0. If X is obtained by applying the filling rules to G and its attached 
radicals (where the radical attached to a is a stopper), then x is a recurrent 
configuration if and only if it has at least one liberty that is a strongly allowed. 

Proof. {=^) Suppose x is recurrent. Let D be the subcluster of G consisting 
of all cells in G whose origin-facing vertex is connected to o by a path of 
vertices with height 3. Also, note that every liberty of x must be attached to 
a vertex in D. Now take S to be the set of all vertices in D and the liberties 
of X (with respect to G). If every liberty of x were weakly allowed, then 
S would form a forbidden subconfiguration. Thus, there must at least one 
liberty that is strongly allowed. 

(<^=) Suppose that f/ is a strongly allowed liberty of x- Apply the burning 
algorithm. Since U is strongly allowed, every vertex in U must burn. Now, 
since U is a liberty, the vertex v to which U is attached burns, and there 
exists a f o-path where each vertex has height 3 — all of those vertices burn 
sequentially. Once the origin burns, its attached radical must burn (since it 
is allowed). The remaining vertices in G burn according to Lemma [9] and the 
proof of Theorem |3| Therefore, x is recurrent. □ 

^Initially, we make a small adaptation in that the radical attached to o can be filled by 
any stopper. 

^Note that the filling rules already require that a liberty be a stopper. 
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5.2 The breakthrough problem 

The filling rules work when each vertex in the graph can only topple once, 
such as in the decorated rooted subtree. However, when the origin o has 
degree d in G, it can topple more than once. We can arrange the toppling 
sequence for the avalanche caused by adding a grain to o so that we only 
topple o if it is the only vertex with more than d grains. At any time we 
are required to topple o, it must have d + 1 grains, and so the unstable 
configuration ip at this time is formed from a stable configuration ip' by 
adding a grain to o. Then, by Lemma [2] and using the arrangement described 
above, we can create a toppling sequence of the form 

0,Vii, . . .,Vinj^,0, V21, . . . ,0, Vkl, ■ ■ ■ ,Vknk 

where for each i from 1 to /c, no vertex is repeated in the sequence f ji, . . . , viui. 

In this way, we can envision an avalanche as a series of "waves" of vertices 
that topple between successive topplings of o. In the case of a tree, no new 
vertices topple after the first wave. This is a consequence of the following 
lemma: 

Lemma 11. Let G he a tree with root a that is a finite subgraph of infinite 
connected d-regular graph T . Let x be a configuration on G such that = 
3. Suppose that v ^ a topples in the avalanche caused by adding a grain to 
X at 0. Let 

O, 1711, ... , Vin^, . . . , O, V21, . . . ,0,Vkl,. . . , Vkn^ 

be the toppling sequence for the avalanche caused by adding a grain to x at 
organized into waves as described above (i.e., a topples if and only if it is 
the only vertex with height greater than d). Then v G {fn, . . . ,'yini}- 

Proof. By induction on the level of v. If v is level 1, then, by Lemma [3} a 
can only topple a second time after v has already toppled; thus, the lemma 
holds. Now suppose the lemma holds for all vertices v whose level is less than 
n, and let f be a vertex of level n. If x(t>) = d and v topples, then it topples 
as soon as its neighbor u on the unique f o-path topples. By the induction 
hypothesis, u topples in the first wave. Therefore, after u topples, we have 
that V has d + 1 grains and can topple, meaning that it is among the vertices 
toppling in the first wave. So the inductive step holds if x(^) = d. However, 
if x(^) < d, then, since v ultimately topples (meaning that it eventually gets 
d + 1 grains), its neighbor u on the unique fo-path must topple twice (see 
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Lemma |4]). By the induction hypothesis, u topples in the first wave. Now, 
consider the configuration ip reached after topphng the vertices in the new 
topphng sequence up to (but not including) the second toppling of o. Note 
that '4){u) < d [u gained at most d—1 grains from toppling of neighbors other 
than V and lost d grains when it toppled) and is the unstable configuration 
reached by adding a grain to stable configuration ip' at o, where ip' is ip minus 
a grain at o. Note that 

O, V21, . . . ,0,ffcl, • • • jVkrik 

is a toppling sequence for the avalanche caused by adding a grain to ip at 
o. By the induction hypothesis, u is among the vertices V21, ■ ■ ■ ,f2n2 that 
topple in the first wave of this new sequence. Also, we apply Lemma |4] to 
show that exactly one neighbor of u topples prior to u in the first wave of 
this new toppling sequence. Thus, when u is scheduled to topple, it has at 
most d—l + l = d grains. But this is impossible. Therefore, u cannot topple 
twice. Thus, f 's toppling implies that xi'^) = which is the case we have 
already handled, and so the induction is complete. □ 

Thus, on trees, any vertex either topples in the first wave of an avalanche 
or never. This is not true, in general, for cacti. To see this, consider the 
following configuration on a subgraph of the expanded cactus. 




1 



Figure 3: A example of a configuration where a vertex (here, a) does not 
topple in the first wave but does ultimately topple. 
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Therefore, the filhng rules described in Section |4] need modification to 
correctly enumerate the recurrent configurations where the cells in a specified 
cluster topple. One possible approach is to divide a finite expanded cactus 
G with origin o into two rooted subtrees Ui and U2 by cutting the edge 
connecting o to its neighbor o' not in the origin cell. Then, a configuration 
X on G can be associated with an ordered pair of configurations (xi, X2) on 
Ui and f/2, respectively, and the process of adding a grain to o and relaxing 
the avalanche that occurs (if there is one) can be performed in the following 
manner: 

1. Let i = 1, X? = Xi, and xl = X2- 

2. Add a grain to Xi*"^"* at o to produce configuration Xi*~^'*'- If Xi'^^'{o) < 
3, set Xi*"* = Xi*"^'*' and xi^ = X2{i — 1), and go to Step 8. 

3. Relax Xi in Gi to a stable configuration Xi*^- 

4. Add a grain to X2* ""^^ at o' to produce configuration X2 ^^X2 ^^'(^^0 — 
3, set X2*'* = X2 ^'''5 and go to Step 8. 

5. Relax X2 to a stable configuration X2^- 

6. Increment i. 

7. Go to step 2. 

8. Stop; the avalanche is finished. The resulting stable configuration is 
associated with the ordered pair (Xi*'',X2^). 

5.3 First-wave topplings and the filling rules 

From the procedure above, it is possible to determine which vertices topple 
in the first wave of an avalanche caused by adding a grain to o. Let x 
be a configuration on finite expanded cactus G with origin o with (xi?X2) 
the ordered pair of configurations on decorated rooted subtrees Ui, U2, as 
defined above. Then the vertices that topple in an avalanche (if any) caused 
by adding a grain to x at o are as follows: 

1. No vertices, if Xi{o) < 3. 
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2. The vertices that would topple in Ui in the avalanche caused by adding 
a grain to Xi o, if = 3 and X2(o') < 3. 

3. The vertices that would topple on xi in the avalanche caused by adding 
a grain to xi at o, plus the vertices that would topple on X2 in the 
avalanche caused by adding a grain to X2 at o', if = X2(o') = 3. 

Therefore, to enumerate the recurrent configurations in which the cells in 
a given cluster C about the origin topple in the first wave when a grain is 
added at o, we must do the following: 

1. If C does not contain the cell attached directly to o, fill Ui according 
to the filling rules, have U2 be a stopper, and ensure that at least one 
liberty is a strongly allowed stopper. 

2. If C contains the cell attached directly to o, fill Ui and U2 according 
to the filling rules, and ensure that at least one liberty is a strongly 
allowed stopper. 

If C is a cluster containing the origin, we say that x satisfies the extended 
filling rules with respect to C if 

1. C does not contain the cell directly attached to o, xi satisfies the filling 
rules with respect to C, and X2 is a stopper, or 

2. C contains the cell directly attached to o, Xi satisfies the filling rules 
with respect to C fl f/i, and X2 satisfies the filling rules with respect to 
Cnf/2. 

We say that such a configuration x satisfies the liberty rule with respect to 
C if at least one liberty in x is strongly allowed. 
We want to prove the following result: 

Theorem 4. Let x be a configuration on G, a finite expanded cactus, and 
let C be a cluster of G. Then x is a recurrent configuration and adding a 
grain to the origin vertex causes exactly the cells in G to topple in the first 
wave if and only if x satisfies the extended filling rules and the liberty rule 
with respect to G. 

We first extend Lemma [TO] to all clusters: 
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Lemma 12. Let G be a finite expanded cactus with origin o, and let C be 
a cluster of G about the origin. If x satisfies the extended filling rules with 
respect to C and its attached radicals, then x a recurrent configuration if 
and only if it has at least one liberty that is a strongly allowed. 

The proof parallels that of Lemma [TOj 

Theorem |4] follows from Theorem [3} Lemma 12 and our characteriza- 
tion of the vertices that topple in the first wave of an avalanche on a finite 
expanded cactus caused by adding a grain at o. 



6 Counting the configurations 

Let {Gi)^i be an increasing graph sequence exhausting infinite expanded 
cactus r, and let e > be given. By Theorem 1, for each positive integer 
n, we can choose a graph Gtv in the sequence such that for every cluster C 
about the origin with n cells, every induced decorated rooted subtree U has 
2 > a;c/ > 2 — e. To examine the limiting behavior of the avalanche dynamics, 
we consider a sufficiently large expanded cactus to not only contain every 
cluster of size up to N but also ensure that every induced decorated rooted 
subtree U satisfies that xjj arbitrarily close to 2. 

Ignoring the liberty rule, the problem of enumerating the total number of 
recurrent configurations on a sufficiently large copy of the expanded cactus 
in which exactly n cells topple in the first wave, divided by Y\d=i ^s{Ui), is 
reduced to finding the total weight of animals of n cells on a rooted 3-regular 
tree where an animals weight is the product of 3 for each internal vertex, 4 
for each medial vertex, and 12 for each terminal vertex (corresponding to the 
number of effective fillings). 

The sum of the weights for rooted animals with n cells is given by the 
term in the power series expansion of the generating functions f{x) satisfying 

f{x) = 12x + 8xf{x) + 3x[f{x)Y 

Also, we can write an equation determining g{x), the generating function 
whose term gives the sum of the weights of all animals with n vertices: 

g{x) = fix) + [fix)]' 

To help determine asymptotic behavior of the coefficients of for fix), 
and gix), we use properties of the generating functions, especially their sin- 
gularities (a good source for the methods involved is [Odl95] ). 
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First, using the quadratic formula, we can solve directly for /(x), as 
follows: 

3x[f{x)f + {8x - l)f{x) + 12a; = 

, 1 - 8x - Vl - 16x - 80x2 

fix) = 7 

ox 

Note that we choose the negative square root to ensure that /(x) is defined 
at X = 0. 

For X > 0, the smallest x for which / is not analytic is since all the 
coefficients of / are nonnegative real numbers, / is analytic for |x| < ^ 
except for ^. |Odl95] 

Because of this, we can use the lemma of Polya: 

Lemma 13 ( |PR87j (see also |HRS75] )). Let /(x) be a generating function 
with the power series expansion 

/(x) = ao + aix + a2X^ + • • • 

with radius of convergence a. Further, suppose that f is analytic for each 
X with \x\ < a except for x = a. Suppose that there exist functions g, k 
analytic and regular in a neighborhood of a, with g{a) ^ 0, such that 

fix) = (l - 9 9ix) + (l - 9 ~* kix) (3) 

where s is not a nonpositive integer, and t = ort < s. Then the coefficients 
{an}'^=i satisfy the asymptotic relation 



Tis) 



The functions / and g can be written in the form of Equation [3} allowing 
for asymptotic expansion. In particular 



fix) 



Vl + 4x 



1 - 20x)~("^) + - (1 - 20x)° 



6x 

Also, 

5x-l n n /16 13 1 



g(x) = „ Vl + 4x(l - 20x)"("5) 



18x2 ' ' V 9 18x 18x2 
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Therefore, the coefficients of the power series expansions of both f{x) and 

3 

g{x) are asymptotic K20"'n~2 for an appropriate constant K in each case. 
Therefore, not taking into account the hberty rule, the fraction of recurrent 
configurations in which topphngs occur in exactly n cells when a grain is 
added to the origin is asymptotic to 

K20''n-I _3 

= K n 2 

5 ■ 20" 

Moreover, we can create a sufficiently large expanded cactus, and the 
results do not change based on which graph sequence was used. This means 
that the cell-wise first-wave critical exponent does not depend on the graph 
sequence we choose. 



7 Accounting for the liberty rule 

Ignoring the liberty rule results in overcounting the number of recurrent 
configurations in which exactly n cells topple by including non-recurrent 
configurations in the total. Fortunately, it is possible to account for the 
liberty rule as well. 

Each configuration x has a given number of liberties. Recall that the 
substitution rules specify that among the liberties, at least one must be 
a strongly allowed stopper. Since ^ of the strongly allowed radicals are 
stoppers as well (see Section |4]), the fraction of the configurations of the 
whole graph ffiting x compliant with the liberty rule is 
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where i{x) is the number of liberties of the cluster configuration x- 
For a cluster C, define 



where x ranges over all configurations of C that satisfy the extended ffiling 
rules, and s maps each function x to the number of terminal cells that have 
a 3-2-2 filling. For each ra, define 0„ by 
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where C runs over all clusters of size n, and i{C),m{C),t{C) denote the 
number of internal, medial, and terminal cells of C . 

Recall that g{x) = Yl'^=o ^nX^, where 6„ is the weighted number of clusters 
about the origin of size n. Define h{x) = Yl'^=o 4'nbnX^\ where (pnbn is the 
number of recurrent configurations in which adding a grain to the origin 
causes exactly n cells to topple in the first wave. 

We claim that 0„ > ^ for each n. To do this, we show that at most ^ of 
the configurations of clusters satisfying the extended filling rules have zero 
liberties. Let % be a configuration of cluster C with zero liberties satisfying 
the extended fiUing rules. We will use the following method to arrive at a 
terminal cell connected by a path of vertices with height 3 to the origin. We 
start at the origin cell. From there, we will construct a connected chain of 
cells, all of which have a vertex that can be connected to the origin by a 
path of vertices with height 3. Given a chain of cells, if the last cell is not 
terminal (satisfying our goal), we choose the next cell as follows: if the last 
cell is medial, then we choose the cell in C farther from the origin; that cell 
is connected to the origin by a path of vertices with height 3 because the 
only valid fillings for the previous cell possible in a zero-liberty configuration 
X are 3-3-1 and 3-3-2, with the 3's facing the two adjoining cells in C. If 
the last cell is internal, then it is either 3-3-3 or 3-3-2, in which case we 
arbitrarily select a cell not already in our chain directly attached to a 3 from 
the previous cell. Because the expanded cactus is treelike and C is finite, 
we must arrive at a terminal cell. This terminal cell must either be 3-2-1 
or 3-2-2. However, we can replace the selected terminal cell with 3-3-1, 3- 
3-2, or 3-3-3, giving the new configuration at least one liberty. For every 
seven zero-liberty configurations on C satisfying the extended filling rules, 
there are five configurations on C with at least one liberty also satisfying the 
extended filling rules. Therefore, at most ^ of the configurations on C have 
zero liberties. Thus, 

' - 12 20 48 

Since this holds for all clusters C of any size, we have ^ Also, it 
follows directly from the definition of that 0n ^ 1 for each n. Thus, 0„ is 
bounded between two positive numbers. 

Therefore, the number of recurrent configurations in which exactly n cells 
topple in the first wave is bounded between ^ and 1 times a function asymp- 

3 

totic to Kn~2. This means that the cell- wise first- wave critical exponent of 
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the expanded cactus is |, and the proof of this paper's theorem is estabhshed. 



8 Open questions and extension to arbitrary 
transitive decorations 

The first open question is to examine what the asymptotic behavior is on 
the expanded cactus for the total masses of avalanches that include more 
than one wave. In decorating a tree, we lose the important property of trees 
that any vertices that topple in an avalanche must topple in the first wave. 
Evaluating how the possibility that new vertices may topple in a second 
or subsequent wave affects the asymptotic behavior of avalanche mass in a 
nontrivial fashion. Thus, for the purposes of gathering evidence for mean- 
field theory, it might be better to consider measures like first-wave critical 
exponent that are easier to determine than the critical exponent for graphs 
like the cactus. 

In Section 1, we noted the following mean- field conjecture 

Conjecture 1. Let T be an infinite d-regular tree and V he a graph quasi- 
isometric to r (including, hut not limited to, a decoration ofT). Then T' has 
the same critical exponent as T, namely |. 

We can propose a similar mean-field conjecture for cell-wise first-wave 
critical exponents: 

Conjecture 2. Let T he an infinite d-regular tree and T' he a graph quasi- 
isometric to r (including, hut not limited to, a decoration ofT). Then V 
has the same cell-wise first-wave critical exponent as the first-wave critical 
exponent ofT, namely |. 

(Recall that, on trees, the first-wave critical exponent is equal to the 



critical exponent by Lemma 11 ) 



Let r be an infinite (i-regular tree, and let F be a connected transitive 
graph on d vertices. (This means that F is fc-regular for some positive integer 
k. We only consider such graphs F for which A; > 1; in particular, we prohibit 
F from being K2.) Then, if we define T' to be the decoration of F where 
each vertex is replaced by a copy of F, we have that F' is transitive and 
k + 1-regular. We can consider the first-wave avalanche cell mass on F'. We 
can then enumerate all possible allowed cell configurations and come up with 
tables similar to Tables [T] and |2] in Section |3] that indicate the following: 



37 



for a given origin cell configuration, what combinations of weakly- 
allowed and strongly-allowed radicals result in a recurrent configura- 
tion 



for a given cell configuration that is at the root of a decorated rooted 
subtree, what combinations of weakly-allowed and strongly-allowed 
radicals result in weakly-allowed and strongly-allowed whole radicals 

We can compute the limiting weak-to-strong ratio x. Using the same methods 
as in the expanded cactus, we can show that for any positive integer n and 
for any cluster C about the origin of n cells, there is a deep enough finite 
subgraph G of F that arc the union of cells, where the ratio of the number 
of recurrent configurations of G to the product of the number of strongly- 
allowed radicals induced by C grows as D" for some D e Z+. This gives the 
exponential factor in the denominator. 

To calculate the exponential factor in the numerator, we use another 
counting method analogous to the extended filling rules for the expanded 
cactus. The same filling rule principles in the expanded cactus apply to the 
general case: each cell in the cluster must be filled in with an allowed cell 
configuration corresponding to its position within the cluster, and apart from 
the liberty rule (which still holds), the requirements of the cluster's induced 
radicals are completely determined by the configuration of the cell to which 
they are directly attached. In the general case, instead of internal, medial, 
and terminal cells, there arc 2^*^^ classes of cells in the cluster, identified 
by subsets of the set {2, 3, ... , d}. For each cell in the cluster, reindex the 
vertices of F so that the closest vertex in the cell to the origin vertex of G 
has index 1. Then the cell belongs to the class 

{i e {2, 3, . . . , d}|a radical is attached to vertex i of the cell} 

We then use logic to calculate for each S C {2,3,...,(i}, the allowed 
cell configurations for a cell in class S and the restrictions on the attached 

radicals, if any. Using these rules, we can compute a weighting factor as 
that gives the number of effective fillings. To count the coefficient in the 
numerator, we use generating functions to count the number of weighted 
clusters. The generating functions f{x) and g{x) in the expanded cactus 
case can be adapted for this purpose. We have 



f{x) = X 



Y: as{f{x)r'-^'\ 

Sc{2,3,...,d} 
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The same argument in the expanded cactus that g{x), the generating function 
for the number of weighted clusters containing the origin, tells us that the 
formula g{x) = f{x) + {f{x)Y applies unchanged to the general case. 

If Conjecture [2] is true, the following conjecture should hold regarding 

Conjecture 3. The generating function g{x) has its closest singularity at the 
origin at x = jy. Moreover, near j^, the dominant term in g{x) is (1 ~ Dx)^ . 

This would immediately imply that the critical exponent of F' is |. 

If Conjecture |3] holds, then the two separate methods of counting configu- 
rations give an exponential factor of D in the denominator and an generating 
function with a singularity at j^. This brings up the following open question: 

Open question 1. Is there a combinatorial argument for why Conjecture 
2 holds in the general case of decoration of a d-regular tree by a connected 
transitive graph with d vertices? 

The two exponential factors in the numerator and the denominator are 
computed using different methods. However, these methods currently rely on 
examining all recurrent cell configurations. The number of cell configurations 
is d^^^, where d is the number of vertices in F, the decorating fc- regular graph. 
Of course, some of these might be forbidden, but exhaustive examination of 
each individual recurrent cell configuration quickly becomes prohibitive for 
larger graphs. Therefore, the best decorating graphs to examine next are 
graphs whose symmetry we can leverage to reduce the number of cases to 
examine, such as complete graphs. 

In the case of a complete graph Kd-, one can see that the allowed cell con- 
figurations are precisely those satisfying that there does not exist a positive 
integer m such that more than m vertices in the cell has height less than or 
equal to m. Moreover, when computing the exponential factor in the numer- 
ator, we can use automorphisms of the vertices in a cell and the attached 
radicals to conclude that, for every m with < m < c? — 1, the weighting fac- 
tors as are equal for every S C. {2, 3, . . . , (i} with 15*1 = m. Therefore, when 
formulating an analogue of the filling rules for these cactus graphs, we can 
classify cells in a cluster based on the number of attached radicals without 
having to reference their positions, just as in the case of the expanded cactus 
(decorated by K-i). 

That said, if Open Question 1 is answered in the affirmative, then we 
can conclude that the broad conjecture that every cactus graph has cell-wise 
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first-wave critical exponent | is true without having to engage in prohibitive 
complex case analysis. 
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